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^ ! 10 at the Slichter mode period. The formahsm makes use of the theory of the 

, ^ 11 Earth's free oscillations based upon an eigenfunction expansion methodology. 

12 We determine the complete response in the form of Green's function obtained 

>: 

OO . 13 from a generalization of Betti's reciprocity theorem. Surficial (surface load, 

^1 14 fluid core pressure), internal (earthquakes, explosions) and external (object 

P\j ■ 15 impact) sources of excitation are investigated to show that the translational 

o: 

CN ■ 16 motion of the inner-core would be best excited by a pressure acting at the 

17 core boundaries at time-scales shorter than the Slichter eigenperiods. 



7 Abstract 

8 We use a normal-mode formalism to compute the response of a spherical, self- 

9 gravitating anelastic PREM-like Earth model to various excitation sources 
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20 1. Introduction 

21 The three free translational oscillations of the inner core, the so-called 



22 Slichter modes ( ISlichterl . Il96ll ). have been a subject of observational contro- 
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versy since the first detection by ISmylid (119921 ) of a triplet of frequencies 

las been sup- 



that he at tributed to t he Shchter modes. This detection 
I 



25 ported by iCourtier et al. 



fcOOOh and 



26 been confirmed by other authors 



Rosat et al 



28 studies (^Crossley 



2006 



Guo et al. 



1992 



Pagiatakis et al 



Hinderer et al 



1995 



2007) but has not 



Jensen et al. 



1995 



20071). Also, it motivated addit i onal theoretica 



Rochester and Peng 



1993 



Rieutord 



2002 



Rogister 



20031 ). The search for the Slichter modes was invigorated by the develop- 



ment of worldwide data recorded 



jy superconducting grayimet ers (SGs) of 



31 the Global Geodynamics Project (IHinderer and Crossley 



3) 



20001). Thanks to 



their long-time stability and low noise level, these relative gravimeters are the 
most suitable instrument s to detect the small s i gnals that w'ould be expe cted 



34 from the Slichter modes (IHinderer et al 



1995 



Rosat et al. 



2003 



2004h 



35 The theory is now bette r understood a nd computation predicts eigen 



periods 
P REM ( 
by 



jetween 4 and 6 h (iRogisterl. 



Dziewonski and Anderson 



20031 ) for the seismological reference 



19811 ) Earth model. A more recent study 



Grinfeld and Wisdom! (|2010| ) states that the period could be much shorter 
because of the kinetics of phase transformations at the inner-core boundary 
(ICB). 

The observation of the Slichter modes is fundamental because, the restor- 
ing force being Archimedean, their periods are directly related to the density 
jump at the ICB. Th is parameter is st i ll poo rly known: by analyzing seismic 
PKiKP/PcP phases, iKoper and Pyld ( l2004l ) fou nd that it should be smaller 
than 450 kg/n a'^, later increased to 520 kg /m^ ( IKoper and Dombrovskayal . 
20051 ). whereas [Masters and Gubbind (|2003[ ) obta ined 820 ± 180 kg/m^ from 



normal modes observation. 



Tkalcic et al. 



( 120091 ) have shown that the un- 



2 



certainties associated with the seismic noise might partiall y explain such 



Gubbins et al 



49 discre pancies for the estimates of the ICB density contrast 

50 (120081 ) have proposed a model with a large overall density jump between the 

51 inner and outer cores of 800 kg/m^ and a sharp density jump of 600 kg/m^ 

52 at the ICB itself. Such a model satisfies both the constraints set by powering 

53 the geodynamo with a reasonable heat flux from the core, and PKP travel- 

54 times and normal mode frequencies. The value of the density jump at ICB 

55 for the PREM model is 600 kg/m^. 

56 This paper aims at evaluating the possible amplitude of the Slichter modes 

57 for various types of excitation sources. 
The s eismic exci tation has be en previously studied by 



CrosslevI (jl992|) and 



Rosat 



SmithI f ll976h . 



(120071 ). They have shown that the best natu- 

60 ral focal mechanism to excite the Slichter mode is a vertical dip-slip source. 

61 The largest magnitude event in the past was the 1960 Chile earthquake with 

62 a magnitude M,„ = 9.6 for the main shock. A foreshock occurred with a mag- 



63 nitude of 9.5 (IKanamori and Cipar 



1974 ). The combination of both events 



64 leads to a seismic source of magnitude My, = 9.8 which would be enough to 

65 excite the Slichter modes to the nanoGal level. However, at such frequen- 

66 cies, the noise levels of SG s are of several nCal, even for the quietest sites 

67 ( iRosat and Hindererl . 120111 ). Earthquakes are therefore not the most suitable 

68 source to excite the Slichter modes to a level sufficient for the SGs to detect 

69 the induced surface gravity effect. 



Surficial pressure flow acting in the core has been considered by 



Greff-Lefftz and Legros 



71 ( 120071 ) as a possible excitation source. In this work, we reconsider the pres- 



sure flow acting in the core using a Green function formalism for a non- 



3 



73 
74 
75 



rotating anelastic PREM Earth model. Then, we investigate the surface 
load and meteoroid impact as possible sources of excitation of the Slichter 
modes. 



76 2. Green function formalism 



77 We consider a spherical non-rotating anelastic Earth model. The dis- 

7B placement u at a point r and time t produced by any body force density f 

79 acting in volume V and surface force density t acting upon surface S can 

80 be written as a convolution of the impulse response G with the entire past 



19981): 



81 history of the forces f and t (iDahlen and Tromp . 

u(r,t)= [ [ G{r,r';t-t')f{r',t')dV'dt' + [ [ G{r,r';t - t')t{r' ,t')d^"^' 

J-oo Jv J-oo Js 



dt' 



82 where r' is the integrated position vecto r. This relation i s infe rred from 

83 Betti's reciprocity relation in seismology ( lAki and Richardsl . Il980l ). Seismic 

84 Green's tensor G of a non-rotating anelastic Earth is given in terms of the 

85 normal-mode complex frequencies Vk = ^kiX + ^jtt) and eigenfunctions by 



G(r,r';t) = J](^z.,)-is,(r)sfc(r')e*"'=*i/(t), 



(2) 



where 3? denotes the real part of the complex expression and H{t) is the 
Heaviside function. 



Tromp and Mitrovical (119991 ) have generalized Betti's reciprocity relation 
to a representation theorem suited for surface-load problems, so that the 
displacement u due to a surface load a located at r' is given by 



u(r,t)= I I (x{r',t')r{r,r';t-t')dJ:'dt' 

J — oo J S 



(3) 



91 where T is the surface-load Green's vector defined by 

r(r, r'; t) = -[G{r, r'; t) ■ V^r') + g(r, r'; t)], (4) 

92 V is the gradient with respect to r', $ is the unperturbed gravitational 

93 potential, and g is 

g(r, r'; t) = ^ J](zz/fc)- Vfc(r)sfc(r')e*'^'^*if(t), (5) 

k 

94 0fc denoting the perturbation of the gravitational potential associated with 

95 the normal mode {sk,(pk}- Green's tensor {G,g} represents the complete 

96 point- source response. 

97 A spheroidal mode of harmonic degree / and order m and radial overtone 
9B number n can be decomposed into three components in spherical coordinates: 

^sT{r) = nUr{r)Yr{eA)r + k-\Vr{r)^9 + k-\Vr{r) ' ^ 



dd ' ' 'sine dcj) 

(6) 



99 where ¥^"^(9,4)) are the re al spherical harmonics of degree / and order 



m 



100 ( iDahlen and Trompl . Il998[ ). k = -^/(mj and r, and are the usual 

101 unit vectors of the spherical coordinates. The associated perturbation of the 

102 gravitational potential has the form 



n^T{^) = nPr{r)Yr{eA). (7) 



103 The eigenfunctions nUp{r), nVi^i^) and nPri''^) functions of the radius 

104 only. Because the model is non-rotating and spherically symmetric, the 21 + 1 

105 eigenfrequencies for each fixed I and n are degenerate into a single eigenfre- 

106 quency that we can therefore denote by The summation over k in Eqs 

107 (|2]) and (|5]) is actually a triple summation over /, m and n. Since the eigen- 

108 frequencies do not depend on m, the summation over m can be performed 



109 usinK the addition theorem for surface spherical harmonics, with the result 



Dahlen and TrompI (119981 ). Eqs 10.28 and 10.34]: 

21 + 1 e*"^'* 



n I 



{nUi{r)Mr')rr'Pio 

+k-'[nUi{r)Mry& -n Vi{r)M{r')&r']Pii 
+k-\Vi{r)nVi{r')^^%^inQ)-^Pa}, 



(8) 



111 where $ = f x and 6 is the angular distance between the receiver at r 

112 and the point source at f' : 



cosO = r • r' = cos cos + sin sin cos (0 — ch'). 



(9) 



113 The Slichter mode is the spheroidal mode of harmonic degree one and 

114 radial overtone number one. For a non-rotating spherical model, the three 

115 Slichter frequencies are degenerate into a single eigenfrequency. As 



Pio(cos O) = cos O 



(10) 



and 

Pii(cos 9) = sin 9, 
the term for which / = 1 and n = 1 in Eq. ([8]) writes: 

iG'i(r,r';t) = — 3? { — [f/(r)f/(r')rr cos 9 



111 



+ ^{U{r)V{r')Y& - V{r)U{r')@v') sin 9 
v2 

+ ]^V{r)V{r')@& cos 9 + ]^V{r)V{r')^^']} (12) 



118 being understood that u = lUi, U 



1^1, V = ^V^ and P = iPi. For PREM, 



the Shchter eigenperiod is 5.42 h f lRogisterl . l2003l ) and the eigenfunctions U, 
V and P are plotted in Fig. lA.ll 

The damping rate depends on the dissipation pro cesses involved. A sum- 
plausible dissip a tion processes is g iven by 



mary o 

mm; 



GuoetaL 



Rosat et al. 



124 cosity has been studied by 



Greff-Lefftz and Legros 



(120071). The role of t he outer core vis - 



Smylie and McMillan 



125 the effect of a mus hy zone at the ICB 



126 magnetic field, by 



Buffett and Goertzl (11995 



127 for the core modes, by 



Crossley et al. 



Pengj ( 



(12000 



) and 



RieutordI (l2002h . 



19971 ). the influence of the 



and the anelastic dissipation 



(I1992I ). Such studies have revealed 



that it is unlikely that the damping factor of the Slichter mode be less than 
2000, corresponding to a damping time of 144 days. In this case, the induced 
surface gravity perturbation should be more easily detectable by SGs. We 
assume, in the following, a quality factor of 2000. 

Using the Green function formalism, we can compute the excitation of 
the Slichter mode by any body or surface forces. 



134 3. Excitation by fluid core pressure 

135 Observational evidence for motions in the core comes from core-sensitive 

136 seismic modes, which have periods smaller than one hour, the free core nu- 

137 tation, which is a rotational mode of nearly-diurnal period, and variations of 

138 the magnetic field that can be related to motions in the core with timescales 

139 larger than one year. Therefore, the dynamics of the fluid core at the Slichter 

140 frequencies lacks observational evidence. Theoretical results suggest that, at 

141 timescales smaller than one day and outside the seismic band, plausible mo- 



7 



tions are to be searched for in the turbulent convection or in the spectrum 
of the core. 

An account of small-scale turbulence driven by convection is given by 



Loperl (|2007| ). The timescale may be less than one day but, because of 



146 the small characteristic length-scales, turbulence is unlikely to excite the 



translation of the w 



lole inner core. 



Valettd (Il989a 



b|) has shown that the inertia-gravity spectrum of an in- 



149 viscid liquid core is continuous and set bounds on it. The bounds de- 



150 pend on both the speed o. 



rotation anc 



Rogister and Valettd ( 120091 ) and iRogisterl ( 120101 ) have suggested that the ro 



squar ed Brunt- Vaisala frequency. 



tational modes might be influenced by the continuous spectrum in which 
they are embedded. In particular, the nearly-diurnal free inner core nutation 
and long-period inner core wobble might be double or even multiple and have 
energy in the liquid core. Pending on the value of the squared Brunt- Vaisala 
frequency in the outer core, the Slichter modes could also be embedded in 
the continuous spectrum. Similarly to what has been found for the two ro- 
tational modes of the inner core, significant motion and pressure variation 
in the liquid core could accompany the Slichter modes. Although a Slichter 
mode with its associated motion in the liquid core should then be considered 
as a single normal mode, we can for simplicity assume that the pressure vari- 



162 ations in the liquid core excite the tra nslation a. 



Buffett 



motions of the inner core. 



( 120101 ) did to investigate the 



This is somewhat the opposite of what 
attenuation of the free inner core nutation: he assumed that the tilt of the 
inner core generates shear layers in the outer core where Ohmic and viscous 
dissipation occur. 
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167 As IGreff-Lefftz and Legroa (120071 ) , we will assume that the pressure at 

168 the CMB takes the following analytical form: 



P'i9, 0, t) = PI 



(13) 



169 where Tq is the starting time of application of the pressure, r is the time 

170 duration of the pressure source and Pq{0, 0) includes three terms of harmonic 

171 degree 1: 



P 



P^Q COS 6 + (P^^ cos + sin 0) sin 9. 



(14) 



172 According to 



Okubo and Endol fll986h and ICreff-Lefftz and LegrosI (120071 ). 



173 the total force exerted at the core boundaries must vanish for the centre of 

174 mass to be kep t fixed. This translates into the so-called Consistency Relation 



175 ( iFarrell 



19721 ) and imposes a relation between the pressures P'^ and P*'^ at 



176 the CMB and ICE, respectively: 



Greff-Lefftz and LegrosI (120071) analytically so. 



elasto-gravitational deformation ( lAlterman et al. 



.2 

c_ pc 
2 ■ 

ic 



ved t 



(15) 



le equations of the 



19591 ) for a simple Earth 



179 model made up of three homogeneous incompressible layers and investigated 

180 the excitation of the Slichter mode by a pressure acting at the outer core 

181 boundaries. In this section, we consider the same simple excitation sources 

182 to test the Green function approach for PREM, which is a more realistic 

183 Earth model. We compute the displacement u by means of Eqs ([1]) and (1T2|1 . 

184 As we are mainly interested in the surface gravity effect, we only need the 

185 radial component of the displacement: 



«,(r,^,0,t) = Agfj{^ 
Ait IV 



[t/(re) / COS e(Pi'o COS e' + (Pi'i COS 0' + Pi'i sin 0') sin 9')dE' 

JCMB 

-Uir-i,) [ cos e(Pi*^ cos 9' + (P{J cos 0' + P^^J sin 0') sin e')dT,']} (16) 

ilCB 

186 At the CMB, c/S' = r^^ sin 9' dO'dcp' and, at the ICB, rfS' = r%smd' dd' d4>' . 

187 Taking Eqs (|9]) and (fT5|) into account, the integration over 9' and 0' gives: 

w,(r,^,0,t) = g?{^SM[f/(re)-f/(r,,)] T e-(*-*')e-[^l'rft' 

(Pi'o COS + Pi'i sin cos + P{=i sin 6 sin 0) } (17) 

188 The integral 

lit) = f e-(*-*')e-[^l'dt' (18) 

J — oo 

189 is calculated in the Appendix. It gives 



m = ^re-('-To),-^^ry.^^ ^ erf(^ + .^)], (19) 

190 where erf denotes the error function. The radial displacement then becomes: 

n.(r,g,0,t) = - J ^(rc)-^(r.e)][^^^-^{/(t)}] 

(P{o COS ^ + Pfi sin ^ cos + Pi\ sin ^ sin 0) , (20) 

191 where 3fJ{ } and 53{ } respectively denote the real and imaginary parts of the 

192 expression between brackets. 

193 The degree-/ gravity variation measured by a gravimeter at the surface of 

194 the Earth is the sum of three terms: the free-air gravity variation owing 

195 to the displacement of the ground in the surrounding unperturbed gravity 

196 field go 

2 

^free = {-^TiGp + -go)U{r,), (21) 
rs 
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197 the inertial acceleration of the ground 



(22) 



and the perturbation of the gravitational attraction 

2 



(23) 



In these expressions, p is the mean density of the Earth. 

The degree- 1 gravity variation measured by a gravimeter is therefore 

^2 



w 1 + 



4Q2 



[PiQ COS 9 + P^^ sin cos + P^^ sin 6 sin (j)] 

- ^{im[~^'U{r,) + ^goU{r,) + ^P(r.)](24) 



We consider a zonal pressure P{g = 150 Pa and compute the induced 
geocentre motion, inner-core translation and surface gravity perturbation for 
both r = 1.5 h and 15 h (Fig. IA.2I) . As the centre of mass is fixed, the geocen- 
tre motion, which is the displacement of the figure cent re with respect to the 



205 centre of mass, corresponds to the surface displacement (iGreff-Lefftz and Legrosl . 



19971). 



For P(q varying from to 1000 Pa and r ranging from 0.1 and 10 h, we 
compute a 2D-map (Fig. IA.3|) of the surface gravity perturbation when the 
exciting pressure has vanished (i.e. t > Tq + r). 

Our results for the PREM model s how excitation amplitudes larg er by 



Greff-Lefftz and Legrosl (l2007f ). who 



70 % than the amplitudes computed by 
used a very simple model made up of three incompressible homogeneous 
layers with a solid inner core, a liquid outer core and a rigid mantle. If 



11 



we consider such a 3-layer model with the average densities of the PREM 
inner core, outer core and mantle, we obtain a period of 3.09 h for the 
Sl ichter mode and excitation arn plitudes in close agreement with the values 
of iGreff-Lefftz and Legrosl (120071 ). They do not provide numerical details for 



the structure of their model but mention that its Slichter period is 3.08 h, 
which is almost equal to the Slichter period of our 3-layer model. Therefore, 
the Green function formalism we have adopted gives the same result as the 
analytical s olution for the degree- 1 deformation of a simple 3-layer model 



222 obtained by 



Greff-Lefftz and Legrosl fl2007f ) 



The difference between the e xcitation amplitudes we obtai n for PREM 



224 and the amplitudes computed by 



Greff-Lefftz and Legrosl (|2007l ) comes from 



the different Earth models. The elasticity of the mantle and inner core, the 
compressi bility of the outer core and the density jump at the ICB all come 



into play. 



Greff-Lefftz and Legrosl (120071 ) mention that the elasticity of the 



mantle perturbs the solutions by 30 %, without any further specifications. 
We have checked that, by making the outer core of our simple 3-layer model 
compressible, with a P-wave velocity of 11083 m/s, the excitation amplitude 
of the Slichter mode is decreased by 40 % with respect to the incompressible 
model. The influence of the compressiblit y and stratificat ion of the core on 



the Slichter mode was also investigated by iRogisterl (l200c 

The perturbation of the surface gravity field is the largest when r is 
smaller than half the Slichter eigenperiod. When r is one fourth of the 
Slichter period, a 10 Pa pressure acting at the CMB, which by Eq. f|T5l) 
imposes a 81 Pa pressure at the ICB, is enough to induce a 10 nGal (0.1 
nm/s^) surface gravity perturbation, which should be detectable by SGs. 



12 



239 4. Excitation by a surface load 

240 A degree—/ surficial mass distribution a* at tlie surface exerts forces 

241 over tlie Eartli in two ways. First, at the interface between the Earth and 

242 the load o"^, the static contact forces give rise to a degree—/ pressure 

P' = 9oct\ (25) 

243 Second, the gravitational attraction of the load a* over the entire Earth is 

244 described by a degree—/ potential 

^ AnG , /(^)'ifr<r, 

<P = 7^7— 7^ < (26) 
+ |^(^)'ifr>r, 

245 Atmospheric pressure models are sampled at 3 h at best. So instead of 

246 using actual data, we write the surface density load in the analytical form: 

a\e, 0, t) = a^,{e, 0)e-[^^ (27) 

247 which is the same as the expression used for the fluid core pressure in Section 

248 131 The degree-one load ctq, like Pq, contains three terms: 

(JoiO, </>) = o"io cos 9 + (cTii cos (p + o"ii sin 0) sin 9. (28) 

249 After integration of Eq. (|3]) over the whole surface, the radial displace- 

250 ment is given by: 

ur{r, 9, 0, t) = [U{r,)go + P(r,)] - 53{/(t)}] 

[cTio COS 9 + ail sin 9 cos cj) + an sin 9 sin (p] , (29) 
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251 and the perturbation of the surface gravity is 



[cTio COS 9 + (Til sin 9 cos + (Xii sin 9 sin 0] 




[-a;2f/(r,) + -(7of/(r.) + -P(r.)] 



(30) 



252 



We use a zonal surface load pressure of 1000 Pa (the surface mass density 



263 is then o"io = Pio/Oo) and we compute the induced geocentre motion, the 

254 inner-core translation and the surface gravity perturbation for two excitation 

255 time-scales (1.5 h and 15 h) (Fig. IA.4I) . When applying a surface load of 

256 1000 Pa during 2r = 3 h, the induced surface gravity perturbation has an 

257 amplitude of 5 nGal (0.05 nm/s^) corresponding to an inner-core translation 

258 of 15 mm and a geocentre motion in the opposite direction with an amplitude 

259 of 0.015 mm. When the excitation time-scale (15 h) is larger than the Slichter 

260 period, the excitation amplitude is smaller. 

261 We also plot the surface gravity perturbation associated to the Slichter 

262 mode excited by a surface load for different excitation time-scales and various 

263 zonal pressure amplitudes in Fig. IA.5I The conclusion is similar to the one 

264 for an internal pressure flow, except that the surface gravity variations are 

265 about 300 times smaller. 

266 In this section, we have estimated the effect of a surface load on the 

267 Slichter mode. Another source that, intuitively, could make the inner core 

268 oscillate is a shock at the surface. Hence, in the next section, we study the 

269 excitation by an object impact on the Earth's surface. 
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270 5. Excitation by an object impact 

271 Stellar objects, such as asteroids or comets, are dragged by the Earth's 

272 atmosphere and reach the Earth's surface at relatively modest velocities, 

273 typically a few tens of km/s. The released energy is co mparable to nuclear 



274 explosions (according to Table 6 of lCoUins et al.l ( 120051 ). from 3.2 to 3.9 10 



275 megatons of TNT, where 1 Mt = 4.2 10^^ J). The collisions between the 

276 Earth and the largest meteoroids, with diameters from hundreds of meters 

277 to several kilometers, blast out the impacting objects, create wide craters, 

278 generally twenty times larger than the diameter of the meteoroids, and melt 

279 terrestrial rocks. Fortunately, such collisions are rare events: statistically, 

280 a 100 to 200 m meteoroid hits the Earth every 1000 years, a 500 to 800 m 

281 meteoroid every 30000 years, and a 5 km meteoroid every 40 millions years. 

282 We believe it is reasonable to assume that both the mass and linear mo- 

283 mentum of the impacting object are negligible with respect to the Earth's 

284 mass and linear momentum, so the orbit of the Earth is not disturbed. Be- 

285 sides, the rotation period and tilt of the rotation axis of the Earth could 

286 be changed by the impact if the angular momentum of the object is large 

287 enough. We, however, consider impactors for which the angular momen- 

288 tum is at least one hundred times smaller than the Earth's; the change of 

289 the Earth's angular momentum is therefore negligible. The major known 

290 meteoroid impacts had such characteristics (Table RTT]) . 

291 Consequently, we reduce the extraterrestrial impact issue to the com- 

292 putation of the equivalent seismic magnitude corresponding to the released 

293 energy. The computati on proposed h e re is b ased on the equations and drastic 



294 simplifications used by iCoUins et al.i (i2005l ). which are summarized below. 
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When an object enters the Earth's atmosphere, it loses its kinetic energy 
296 through deceleration and a blation. The r ate o f change of th e velocity v is 



297 given by the drag equation (jColhns et al 



2005 



Melosh 



19891): 



dv 



SpzCo 2 

-V , 



(31) 



298 where z is the altitude, Co is the drag coefficient, taken equal to 2, and pi 

299 and Lq are, respectively, the impactor density and diameter. By assuming 

300 an exponential atmosphere. 



p{z) = poe-^/^. 



(32) 



where H = 8 km is the scale height and po =1 kg/m^, the velocity of the 
impactor as a function of altitude is given by: 

3p{z)CdH' 



v{z) = Vq exp 



(33) 



4pjLo sin a 

where a is the entry angle and vo, the velocity at the top of the atmosphere. 
On its trajectory down to the ground, the impactor goes through t he in- 
creasing atmospheric pressure and, possibly, breaks up. iColhns et al.l (|2005[ ) 
have established an empirical strength-density relation to estimate the yield 
strength Yi (in Pa) 

logio Yi = 2.107 + 0.0624v^ (34) 
and give an approximate expression for the altitude of breakup z*: 

Y 



-H[\n{ 



^Povt 



+ 1.308 - 0.314// - 1.303v/l -//], 



where 



4.07 



CpHY, 
PiL^vf sin a ' 
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(35) 



(36) 



310 and Vi is the impactor velocity at the surface. 

311 Equation f p5|) holds provided that If < 1. Otherwise, and more rarely, 

312 the object does not break up and the velocity at the impact is given by 

313 Equation fl33]K 



314 Following 



315 ( IChyba et al. 



Collins et al 



1993 



Melosh 



( 20051). we use the approximative pancake model 



198ll ) to describe the disintegration of the me- 



316 teoroid. Let us denote by Zb the airburst altitude, which is the altitude of 

317 complete dispersion of the fragments. According to the simplifying assump- 

318 tions of the pancake model, Zb is given by 



Zb = z"- -2Hln{l + 
where i is the dispersion length scale: 

i = Lq sin a 

and t 



2H 



(37) 



(38) 



le pancake factor /„ is between 2 and 10. We shall adopt 



Collins et al 



(120051 ) 's value of 7. If Zb > 0, the airburst occurs in the atmosphere, there is 



neither impact nor associated seismic event. If Zb < 0, the fragments are not 
dispersed when they collide with the ground and the impact velocity is: 



t;(2: )exp{(-- 



4 piLl sin a 2>P' 



(32 + (^I'e^*/^ 



The remaining kinetic energy at the moment of impact is 

TT 
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(39) 



(40) 



After the imp act, a fraction fc., o i E^ i s radiated as seismic waves. Exper- 
imental data ( ISchultz and Gaultl . 119751 ) provide ks G [10~^,10~^]. We will 
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327 take kg = 10 ^. The seismic moment being given by 

Mo = 2-^ksE,r, (41) 

328 where the stress release Acr^ ^ 3 MPa and the rigidity /x = 30 GPa, the 

329 seismic magnitude is then: 

M^ = ^logio(Mo)- 10.73 (42) 

330 with Mq in dyn.cm (1 dyn.cm=10~'' N/m). 

331 We compute the magnitude My^ for the different meteoroid impacts of 

332 Table lA.ll and, in the same table, report the surface gravity perturbation 

333 associated with the translational excitation of the inner core. Note that the 

334 source is represented in terms of moment tensor by three orthogonal force 

335 couples (spherically symmetric explosion) and not as a vertical force. Indeed, 

336 the shock pressure would reach hundreds of gigapascals and the impact en- 

337 ergy would vaporize the rocks and cause a spherically symmetric explosion, 

338 as observed from the spherical shape of known craters. 

339 For the biggest meteoroid, the surface excitation amplitude of the Slichter 

340 mode is 0.0067 nm/s^, which is less than the detection threshold of 1 nGal 

341 (=0.01 nm/s^). To determine what kind of impact would be necessary to 

342 excite sufficiently the Slichter mode so that it is detectable in surface gravity 

343 data, we compute the magnitude and gravity perturbation Ag for various 

344 ranges of density pi, velocity Vq, diameter Lq of the object and for seismic 

345 efficiency kg varying between 10~^ and 10^^. The resulting maps are plotted 

346 in Fig JA.6[ The shaded areas correspond to My, larger than 9.7, which is the 

347 magnitude required for the surface gravity effect to reach the nGal detection 

348 threshold for a surficial explosive moment source (Fig. IA.7p . 
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349 We see from Fig JA.6l that to produce a seismic event of magnitude larger 

350 than 9.7, the size, density or velocity of the impacting object should have 

351 unrealistic huge values. However, the value of the seismic efficiency ks has 

352 a direct impact on the equivalent seismic magnitude. For instance, with a 

353 seismic efficiency of 10^^ instead of 10"'*, a meteoroid similar to the one 

354 which produced the Chicxulub crater in Mexico would be able to induce a 

355 seismic event of such a magnitude. Of course, the consequences would have 

356 been devastating. 

357 We conclude that the surficial seismic events, including extra-terrestrial 

358 object impacts and explosions, are not efficient to make the inner-core oscil- 

359 late at the Slichter frequency with an amplitude large enough to be observed 

360 at the surface. The reason is the same as for earthquakes, i.e. the excitation 

361 amplitude is directly linked to the seismic magnitude and the radial eigen- 

362 functions of the Slichter mode are constant and close to zero in the mantle 



(Crosslev. 


1992; 


Roeister. 


2003; 


Rosat, 


2007) 



364 6. Conclusions and perspectives 

365 We have investigated the excitation of the translational free motion of the 

366 inner core by a pressure due to a flow in the outer core and acting at both the 

367 ICB and CMB, by a surface load, which can be associated to atmospheric or 

368 oceanic loading for instance, and by the collision between the Earth and a 

369 stellar object. Our conclusion is that the Slichter mode would be best excited 

370 by a pressure acting at the core boundaries at time-scales shorter than half 

371 the Slichter eigenperiod. 

372 For the pressure source at the ICB and CMB and the loading source at the 
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surface, we have considered Gaussian functions of time. More complicated 
sources should be considered, in particular stochastic forces produced by 
some turbulent flow in the core or at the surface. The stochastic excitation, 
be it oceanic and atmospher ic , of norma l modes has been studied for instanc e 



by 



Tanimoto and Um 



(119991 ) 



Tanimotol (11999 



20071 ) and lWebbI (l2007l . l2008h 



However, the time scale for the Slichter mode is larger than for the other 
seismic normal modes, whose eigenperiod is shorter than 1 hour, so we should 
consider a theory different from the Kolmogorov theory of turbulence. 

We have considered an analytic expression for the surface pressure as a 
source for the excitation of the Slichter mode. A more realistic approach 
should be based on actual atmospheric data from space correlation of world- 
wide barometers or from weather institutes (ECMWF, NCEP...), provided 
the data are available at a time resolution higher than the Slichter period. 
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386 Appendix A. Evaluation of the time integral I(t) (Eq. I18|) 

387 We evaluate the following integral: 



CO 



388 We put b = and introduce the variable x = I{t) becomes: 

/b pb 
-oo J — oo 

389 We perform the change of variable y = x + i^: 

I(t) = re'''^' ' e-(^-'^)(^+*^)c/l/ 

J — oo 

2 2 /•^+«^ 



390 The integral j^'^ e dy is the complementary error function erfc(z) = 

391 1 — erf (2;), where 

erf (2) = ^ / e-y^dy. (A.l) 
V^r Jo 

392 So the time integral I{t) is given by: 



I(t) = ^^e*K*-To)e-.VV4[i + erf(^-^ + z^)] 
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Table A.l: Some meteoroid impacts on the Earth continental crust and oceanic crust. The 
impact angle is supposed to be 45 degrees and the impact velocity is 20 km/s. 





Date 


Diameter 


Density 






Location 


(AD or My BP) 


(m) 


(ker/m^) 




(nm/s^) 


± Ullg, LlolwX J7 11 C Udll 


ions AD 




9700 


No impact 


kjl lyCl la 






I 1 U JV \ 






Ries Crater 


15.1 ± 0.1 


1500 


2700 


7.4 


3.9 10-6 


Germany 






(rock) 






Rochechouart 


214 ±8 


1500 


3350 


7.5 


4.9 10-6 


Prance 






(stony-iron) 






Chesapeake Bay 


35.5 ±0.3 


2300 


2700 


7.8 


1.4 10-5 


USA 






(rock) 






Chicxulub 


64.98 ± 0.05 


17500 


2700 


9.6 


6.7 10-3 


Mexico 






(rock) 
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Figure A.l: Eigenfunctions of the Slichter mode iSi for the PREM model. The vertical 
axis is the radius normalized by the Earth's surface radius r^. U and V are the radial 
dependence of the displacement given by Eq. ([5]) . P is the perturbation of the gravitational 
potential. The normalization of the eigenfunctions is such that U{rs) — I la. 
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Figure A. 2: Effects of the excitation of the Shchter mode by a fluid pressure acting at 
the CMB for two different excitation time-scales: r = 1.5 h (solid line) and r = 15 h 



(dashed line). This figure is similar to Fig. 3 of iGreff-Lefftz and Legrod (j2007h but we 
have applied a Green tensor formalism to the PREM model, (a) Degree-one zonal pressure; 
(b) geocentre motion; (c) inner-core translation; (d) surface gravity perturbation. 
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Figure A. 3: Surface gravity perturbation induced by the Slichter mode excited by a fluid 
pressure acting at the CMB for different excitation time-scales and various zonal pressure 
amplitudes. The vertical dotted lines correspond to one fourth of the Slichter period, one 
half of the Slichter period and the Slichter period. 
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Figure A. 4: Resulting effects of tlie exeitation of the Slichter mode by a zonal surface 
load a{t) ~ (Toe~'~i^l for two different excitation time-scales: r = 1.5 h (solid line) and 
r = 15 h (dashed line), (a) Degree-one zonal pressure effect; (b) geocentre motion; (c) 
inner-core translation; (d) surface gravity perturbation. 
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Figure A. 5: Surface gravity perturbation induced by the Slichter mode excited by a surface 
load for different excitation time-scales and various zonal pressure amplitudes. 
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Figure A. 6: Seismic magnitude as a function of the impactor parameters and seismic 
efficiency kg. The shaded area corresponds to moment magnitudes larger than 9.7, i.e. to 
induced surface gravity perturbation larger than 1 nGal. 
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Figure A. 7: Surface gravity perturbation induced by the Slichter mode as a function of 
the moment magnitude of a superficial energy release (explosion or object impact). 
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